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We present experimental data on ultracold neutron (UCN) production in solid oxygen from 5 K 
to 50 K. The experiment used 3 cells of different dimensions in order to study the excess UCN 
elastic scattering in S-O2 due to source density non-uniformities. We present a combined analysis 
with the aid of Monte-Carlo simulations which contain the detailed physics of UCN transport as 
well as cold neutron attenuation. The results validate the UCN production and upscattering cross- 
sections extracted from the dynamic structure function derived from inelastic neutron scattering, 
and illustrate the role of magnetic scattering in enhancing UCN production in S-O2. 
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As superthermal ultracold neutron (UCN) sources are 
currently under construction around the world, using ei- 
ther deuterium or superfluid helium as converter mate- 
rials, the investigation of other materials is a potentially 
interesting topic to a broad scientific community. In par- 
ticular, oxygen, from a theoretical point of view, seems to 
be a promising alternative, as it has a lower nuclear ab- 
sorption (compared to deuterium), is commercially avail- 
able with high purities, and has sufficiently low Fermi- 
potential, although its UCN production cross section is 
slightly smaller than that of deuterium. Detailed study 
of this new material could yield new scientific insight. 
Unfortunately, experimental data presented in the past 
few years [l|-|3| are inconclusive and even contradictory. 
Despite the fact that these experiments were carefully 
designed to provide direct comparison of UCN produc- 
tion from the materials under study, interpretation of the 
results was difficult, mainly because of the intrinsically 
different thermal properties of these cryogenic solids that 
limit the knowledge of UCN extraction. 

In a typical UCN production experiment, the observed 
UCN count rates are a convolution of (a) UCN produc- 
tion, (b) UCN loss during extraction from the source, (c) 
UCN transport from the source to the detector, and (d) 
the variation of the detector efficiency due to the differ- 
ent UCN spectrum as a result of the velocity boost that 
all UCN acquire upon exiting the source material into 
a vacuum guide tube. The UCN transport can be sim- 
ulated using detailed Monte-Carlo tools, given that the 
physics of UCN transport in a vacuum guide is known 
relatively well. The loss of UCN in the process of extrac- 
tion from the production source, on the other hand, is 
highly dependent on the physical condition of the source. 
Any variations in material density due to cracks or voids 
can significantly change the UCN mean free path (mfp). 
Even though visual inspections were frequently made 



during experiments, most assessments of the solids grown 
in-situ are qualitative at best and might not be directly 
applicable to UCN, due to the fact that the wavelength 
of UCN is an order of magnitude shorter than that of 
optical photons. 

Furthermore, previous analyses often lack sufficient 
treatment of the evolution of the incident cold neutron 
(CN) flux to reliably extract the physics of UCN pro- 
duction in the source material under study. To better 
understand the UCN production mechanism, we have de- 
coupled the CN production from the UCN production by 
placing the experiment on a CN beam line, even though 
the experiment takes a significant cut in the available CN 
flux. The UCN production source is often made thick in 
order to increase the probability of downscattcring a CN 
into a UCN. In such a source, we expect that CN ex- 
perience strong elastic Bragg scatterings that are highly 
energy dependent. Unfortunately, this piece of physics 
is often neglected, most likely because the standard neu- 
tron transport codes (e.g., MCNP Q) do not contain 
readily available data files that include coherent scatter- 
ing [5|. There remain many difficulties in ascertaining 
coherent understanding among many UCN production 
experiments. In this letter, we present an analysis of 
UCN production data in solid oxygen (S-O2), that at- 
tempts to take into account of all the effects mentioned 
above. Experimental data were taken using three cells 
of different dimensions to give us a handle on the mfp of 
UCN. Combined with additional data on inelastic neu- 
tron scattering, we demonstrate how the analysis sheds 
light on the mechanism of UCN production in S-O2 in 
low temperature phases at the saturated vapor pressure. 

The experimental data presented here were taken at 
Lujan Center neutron flight path 12 at LANSCE during 
the summer of 2008. Details of the experimental setup 
was reported in The performance of the apparatus 
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was bench-marked using solid ortho-deuterium in the tar- 
get cell maintained between 5 K and 22 K. Because of the 
low saturated vapor pressure of S-O2, vapor deposition 
techniques used to attain optically transparent solid deu- 
terium do not work for oxygen. Instead, gaseous oxy- 
gen (at 99.999% purity) was first condensed into liquid 
around 60 K in the target cell, and then slowly cooled into 
solid phases. Unlike deuterium, oxygen has three distinct 
solid phases at saturated vapor pressure: 7 (44 ^ 55 K, 
cubic lattice, paramagnetic), /? (24 ^ 44 K, rhombohe- 
dral lattice, short range anti-ferromagnetic(AF)), and a 
phase ( < 24 K, monoclinic lattice, 2-D AF). The strong 
magnetic interaction between O2 molecules (each with a 
non-zero electronic spin S = 1) plays a significant role 
in the abrupt changes of lattice structures among these 
low temperature phases. We use Bragg edges in the to- 
tal cross section (obtained from CN transmission data) 
to verify these lattice structures. 

UCN production was studied using three target cells of 
different length (small: 1.14 cm, medium: 3.56 cm, large: 
8.61 cm), all with the same inner diameter of 6.40 cm. 
The background subtracted UCN production signals, S, 
(the UCN count rate normalized to the incident CN cur- 
rent as defined in [^) from the three cells are plotted 
in Fig. 131 The target cell was cooled slowly from liq- 
uid at 60 K, through 7, /3, and into a-02 at 5 K over 
the duration of 5 days (following the recipe in Q). The 
major challenge arises from the reduction of the molar 
volume of S-O2 by 12.5% upon cooling through the three 
solid phases at saturated vapor presure. A discrete 5% 
change at the 7-/3 transition at 44 K is the bottle- 
neck to attaining large-sized S-O2 cryocrystals in the /3 
and a phases. Even with our best efforts in temperature 
control, subsequent faster warm-up and cool-down cycles 
(over 10 hours) did not show any noticeable changes in 
the UCN production. At the first glance, the data pre- 
sented in Fig. |3] beg for explanations for (a) their non- 
superthermal behavior, (b) their non-linear volume de- 
pendence, and (c) the peculiar production peak at the 
(5 — a transition. The remainder of the paper is directed 
to explaining these atypical behaviors. 

Due to the distinct magnetic orderings, we expect the 
magnon excitations to vary among the three solid phases, 
even without detailed a priori understanding of the un- 
derlying physics. A follow-up measurement [9| of the dy- 
namic structure function, 5(|(5|,aj), of S-O2 was carried 
out using the Disk Chopper Spectrometer at the NIST 
neutron scattering center in 2009. As described in 
we can extract the differential cross section for UCN pro- 
duction and the upscattering using S{Q, lu), regardless of 
the properties of the modes of excitation, i.e.. 
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where uj = E — Eucn is the energy transfer, and to„ is 
the mass of a free neutron. Integrating Eq. [l]over the in- 
cident CN flux, we then get the spectrum-averaged UCN 
production rate; integrating Eq. [2] over all final phase 
space, we get the UCN upscattering cross section. Tables 
of S{Q,ui) were measured at 28 temperature points be- 
tween 5 K and 47 K. The extracted UCN production rate 
and the upscattering cross section are plotted in Fig. [TJ 

The anomalously high production of UCN in the /3 
phase can be explained by its large UCN production rate 
(as shown in Fig. [!}, which primarily originates from 
exciting the soft diffuse mode clustered around a fixed 
Q = 1.3 A and extending from cj = ~' 8 mcV 0]. 
These modes are associated with the short range AF 
clusters present in the /3-02 and should be related to 
the interesting cross-over behavior of the geometrically 
frustrated spins on a triangular lattice before the spins 
arc fully ordered. The physics describing these excita- 
tions has not yet been fully formulated. Nevertheless, 
the inelastic scattering data provides the handle we need 
to proceed to understand the results of our UCN produc- 
tion experiment. 

In a canonical superthermal source, where the loss 
mechanism is dominated by the upscattering loss, the 
ultimate achievable UCN density can be calculated from 

Pucn = PXT= {NadownCkcN) X ^ = I* 

is the ratio of the downscattering to the upscattering 
cross sections that gives rise to the famous superther- 
mal temperature dependence oc exp(_E/fcT) (see the in- 
set in Fig. [1} . The superthermal behavior persists over a 
wide range of temperatures even when the mechanism 
of scattering changes abruptly between phases. Note 
that this behavior is only evident with a large enough 
source, in which the UCN mfp is limited by the upscat- 
tering. With the reduction of upscattering by deep cool- 
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FIG. 1: Temperature dependent UCN production rate and 
UCN upscattering cross section in S-O2. Inset: The ultimate 
achievable UCN density by taking adomn/cTup- 
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ing, other sources of finite loss would become comparable 
in strength and finally limit the achievable gain. 

Next, the UCN transport in the UCN guide system can 
be understood using Monte-Carlo simulations @. How- 
ever, a UCN shutter was installed in the middle of the 
experiment leading to a reduction of the overall trans- 
port efficiency for the data set taken with small and large 
cells. Lacking sufficient data to isolate this effect, we 
introduce a free parameter, f guide, to rescale the UCN 
signals for the small and large cells uniformly up, in or- 
der to compare with data taken using the medium sized 
cell. The UCN spectrum produced inside the source is 
described by the v'^dv distribution, with the velocity v ex- 
tending beyond the critical velocity of UCN. The UCN 
guide serves as a spectrum filter, allowing the low ve- 
locity UCN (with energies below the Fermi potential of 
the guide, after gaining an energy boost of 87 ncV ex- 
iting the oxygen source and entering the guide vacuum) 
to be transported with f 00% efficiency and filtering out 
the high velocity neutrons (very cold neutrons, VCN). 
However, because of the large angular acceptance and 
short length of our UCN guides, there remains a signifi- 
cant fraction of VCN reaching the detector. A reasonable 
neutron source distribution (derived from a Monte-Carlo 
simulation of UCN transport through the guide system) 
is presented in Fig. [5^ to account for these transport ef- 
fects in our apparatus. 

Finally, the production of UCN inside the finite sized 
cell should not be spatially uniform because the CN fiux 
is attenuated as it passes through the cell, mainly by 
elastic scattering. In our initial analysis, we scale the 
production rate, Fig.([T]), by an exponential decay func- 
tion along the length of the cell, leaving the CN mfp as a 
free parameter to be determined by the best experimen- 
tal fit. Note that this treatment assumes that the CN 
attenuates at the same rate independent of the energy. 
The model was later refined to include the evolution of 
the energy spectrum due to the strong Bragg scattering 
that is highly dependent on the CN wavelength. 

We build a model with a UCN velocity distribution ac- 
cording to Fig. [2^, and the volume distribution according 
to the CN attenuation, using the temperature-dependent 
UCN production rate (Fig. [T]). We then allow the indi- 
vidual source UCN to propagate uniformly through An 
with the physics of UCN diffusion described by an excess 
elastic mfp to account for the material non-uniformity 
(Otherwise, the UCN elastic mfp is infinite in a uniform 
O2 solid). The process of UCN diffusion also includes 
the energy-dependent upscattering loss and the nuclear 
absorption loss. The temperature dependence of den- 
sity variation is included in the model. In addition to 
this excess elastic mfp, we add in additional 0.5 b and 
2.5 b to the UCN elastic scattering in /3-02 and 7-O2 
respectively, as suggested by the increased total cross- 
section measured with long-wavelength neutrons below 
the lowest Bragg cutoff. This additional scattering prob- 
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FIG. 2: (a) The UCN velocity distribution inside the source 
used in the Monte-Carlo simulation, (b) Reduced a-ud 
guide transmission rescale factor f guide to correct for the re- 
duced transmission due to the installation of a UCN shutter, 
which affects the data sets collected using the small and large 
cells. Contour plots of the reduced with varying upscat- 
tering scale factor, fup, and the CN mfp (c), and the UCN 
excess cross section (d). 



ably originates from paramagnetic scattering. The model 
assumes that the excess elastic scattering (due to the 
physical cracks) is fixed among all phases. It is probably 
well justified, considering that the change of the lattice 
structure and the molar volume from /? to a is minimal, 
and care were taken to avoid thermal shocks during the 
experiment. The CN attenuation is also treated the same 
in all phases, as the difference of the total cross section in 
a and /? is minuscule, due to their similar lattice struc- 
ture and lattice constants (but it is not quite adequate 
to extend this assumption to the 7 phase). 

We then perform Monte-Carlo simulations using this 
basic model to calculate UCN yields from the 3 different 
cells, scanning through the full parameter space allowing 
variations of f guide, CN mfp, and the upscattering nor- 
malization factor fup that translates the upscattering loss 
(in Fig. [T|) into the absolute cross section. The reduced 
for each set of parameters was obtained by allowing 
the overall normalization, TV, to be adjusted until the 
combined of the data from 3 cells. 
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is minimized for each temperature point independently 
(as N contains no physical significance). The cell index 
runs from small, medium, to large. Results of this study 
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Temperature (K) FIG. 4: Simulated results of the UCN yield from S-O2 ex- 

tracted from (a) downstream of the CN flux (front extrac- 
tion), (b) perpendicular to the CN flux (side extraction), and 
FIG. 3: Best combined flt of the UCN yields from solid O2 in (g) upstream of the CN flux (back extraction), 
cells of three different lengths. 



arc shown in Fig.[2j The best fit indicates that the trans- 
mission was reduced to about 62% with the installation 
of an open UCN shutter, the upscattering scale factor is 
around 1.9, the CN mfp is about 14 cm, and the excess 
UCN elastic cross section is below 1 barn. The best fit 
with these optimum parameters to the experimental data 
is shown in Fig. [3l Note that not only the temperature 
dependence, but also the non-linear volume dependence 
of the UCN yield are reproduced in the simulation. 

The same simulation predicts different UCN yields for 
UCN extracted from different surfaces of our cylindrical 
cell (see Fig. |4]). In our experiment, the CN flux enters 
the cell horizontally from one surface and exits from the 
opposite surface. In what we shall term "front extrac- 
tion", the UCN produced are extracted downstream of 
the CN flux. The UCN acceptance into the guide sys- 
tem is largest for UCN born in the volume slice close to 
the downstream window. On the other hand, the CN are 
significantly attenuated downstream, leading to a smaller 
amount of UCN produced in this volume. This is why 
the large cell does not produce the largest UCN yield 
in our front extraction configuration. The size of this 
effective volume, in which the UCN can be efficiently ex- 
tracted, scales with the total mfp, combining effects of 
elastic scattering, upscattering and nuclear absorption. 
The expected superthermal temperature dependence is 
only evident when the UCN mfp is limited by the up- 
scattering, and is smaller or comparable to the small- 
est dimension of the cell. In a-02, both the production 
and upscattering are smaller, and thus the upscattering 
mfp is already much larger than the cell length, result- 
ing in a very weak temperature dependence as shown 
in both the experimental data and the simulations. On 




FIG. 5: (a) Total cross section of CN scattering in S-O2. (b) 
Simulated results of the CN flux attenuation (along the depth 
of the cell Z) and the energy spectrum evolution, (c) UCN 
production rate weighted by the modified CN spectrum, (d) 
Linear fits to the reduced UCN production rate. 

the other hand, the number of UCN extracted from the 
upstream window (back extraction), and from the side 
(side extraction, including vertical extraction) increases 
monotonically with the increasing volume of S-O2. Most 
noticeably, the expected temperature dependence from a 
superthermal source is restored in a large cell with UCN 
extracted from the side of the cell (into a UCN guide 
oriented perpendicular to the CN beam direction). This 
explains the seeming contradiction between the experi- 
mental data reported in Q and in this paper, as a purely 
geometrical effect. 
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As the CN flux is attenuated in the solid, its energy 
spectrum also evolves as a result of the energy-dependent 
cross-section (see Fig. [5^). This modifies the spectrum- 
averaged UCN production rate used in the simulation 
discussed above. To study this effect, we developed an- 
other Monte-Carlo simulation (based on Geant4) for CN 
tracking, in which the CN elastically scatter isotropically 
with probabilities determined by the energy dependent 
total cross section, which is in fact dominantly elastic. 
Since the production of UCN in the source does not de- 
pend on the direction of the CN flux (incident or scat- 
tered), we record the usable flux at each point in space 
without taking the projection along any specific normal 
vector. The resulting flux attenuation over the depth of 
the cell, Z. and spectrum evolution are shown in Fig.[5]3. 
With this information, we proceed to calculate the spec- 
trum averaged UCN production rate at each volume slice 
along Z, for individual temperature points of interest 
(Fig. [5t) . The reduction of the UCN production rate over 
the length of the cell can be fitted by a linear function 
(the first terms of the Taylor expansion of the exponen- 
tial function). If we were to recast the reduction of the 
production rate as solely due to the CN flux attenuation, 
we arrive at a mfp of 12.5 cm, which is slightly smaller 
than the best value of 14 cm (obtained in the previous 
analysis where the CN mfp is left as a free parameter). 
This removes one free parameter required to perform the 
combined analysis described earlier. The revised anal- 
ysis yields a larger x^, leaving the best values of f guide, 
fup, and the excess UCN elastic cross section the same. 
However, the constraint on the UCN elastic cross sec- 



tion is not tight enough to reach any definite conclusion 
about the limiting size of a practical UCN source made 
of S-O2. The success of the fit validates the UCN produc- 
tion cross section extracted from the dynamic structure 
function. We conclude that the UCN production in S-O2 
is highly enhanced by inelastic magnetic scattering, in 
particular, the large UCN production observed in the /3 
phase can be attributed to the excitations of a geomet- 
rically frustrated spin system. This work was supported 
by NSF grant 0457219, 0758018. 
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